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Summary 

Analytic representations are developed for the discrete 
blade deflection and the continuous tip static pressure 
fields in a stationary reference frame. Considered are the 
sampling rates equal to the rotational frequency, equal to 
the blade passing frequency, and for the pressure, equal 
to a multiple of the blade passing frequency. For the last 
two rates the expressions for determining the nodal diam- 
eters from the spectra are included. A procedure is 
presented for transforming the complete unsteady 
pressure field into a rotating frame of reference. The 
determination of the true flutter frequency by using two 
sensors is described. To illustrate their use, the developed 
procedures are used to interpret selected experimental 
results. 


Introduction 

Over the past few years several experimental programs 
have been completed that used stationary transducers to 
measure the flutter vibrations of rotor blades. These 
experiments were conducted on full-scale engines 
installed in an altitude test facility. Two categories of 
stationary transducers were used: specially designed 
optical displacement sensors, and commercially available 
pressure transducers that could be flushmounted in the 
engine casing or incorporated in a traversing probe. 

The optical blade displacement measurements are 
described in detail in reference 1 . Briefly, this method is 
based on the measurement of reflected light from the 
blade tips. The incident light from a light source and the 
reflected light from the blades are transmitted by means 
of a fiber optics bundle and a lens system. The resulting 
instantaneous blade tip positions are recorded on a 
magnetic tape in the form of a pulse train. In addition, a 
reference pulse (such as a once-per-revolution signal from 
the rotor shaft) must be recorded to provide the base 
relative to which the instantaneous blade positions can be 
measured. Since pressure measurement using miniature 
high-response transducers is fairly common, it will not be 
described in this report. 

These measurements are attractive because all blades 
can be viewed with a single sensor. In contrast, for blade- 
mounted sensors, the number of sensors must equal the 
number of blades in order to measure the response of all 
of the blades. In addition, with stationary sensors the 
complexities and the limitations associated with the signal 
transmission from a rotating component are completely 
avoided. On the other hand, the frequency content of the 
signals is much higher with these sensors, and the quan- 
titative analyses of these measurements are necessarily 
more complicated. For the displacement data the high- 
frequency requirement is necessary in order to be able to 


reproduce and locate precisely the blade pulses. 
Similarly, for the pressure data this requirement is 
necessary in order to be able to reconstruct the steep 
pressure variation (both, the unsteady and the steady 
state) when crossing from the pressure to the suction 
surface of a blade. These requirements were met by 
recording the signals on frequency-modulated tape at a 
high speed and later reproducing the signals during 
digitization at a low speed. Typically, the tape speed 
reduction ratio was 64 and the resulting effective rate of 
digitization was of the order of 1 MHz. Subsequent data 
processing and analyses were performed on a mini- 
computer equipped with extensive storage peripherals. 

Stationary sensors were used to measure displacements 
and pressures during flutter in several previous works 
(refs. 2 to 5). Special procedures had to be developed to 
analyze these data. This report documents the analysis 
methods developed in the past by using mathematical 
derivations and presents a new procedure that allows the 
complete pressure spectrum to be interpreted and the 
complete unsteady pressure field to be transformed into a 
rotating frame of reference. The derived analytic expres- 
sions are particularly essential for the latter procedure. 
Since the problem is by nature discrete, because there are 
a discrete number of blades on the rotor, all of the 
derivations are performed in the discrete domain. From 
another point of view, because of this property, the 
digital spectral analysis methods are ideally suited for 
analyzing experimental data. During the derivations the 
only restriction imposed on flutter frequency is that it be 
less than the blade passing frequency. 

The experimental results presented in the report were 
all obtained for shrouded fan rotors operating at part 
speed and at high incidence. The inlet pressure and 
temperature were raised significantly above the atmos- 
pheric condition and the relative tip Mach number was in 
the transonic range. 

During the course of preparing the report several 
improvements in the mathematical derivations were made 
as a result of discussions with Donald C. Braun of the 
NASA Lewis Research Center. 


Symbols 

A amplitude associated with a blade (real) 
a frequency variable in engine orders (eqs. (18), 
(59), (90)) 

B constant (eq. (3)) 

C aerodynamic coefficient 

c blade chord 

d displacement in plane of disk (fig. 7) 

E engine order 

F defined by eq. (44) (complex) 


1 



/ pressure amplitude (real) 

/ defined by eq. (43) (complex) 

G cross-power amplitude (complex) 

i imaginary unit 

j integer, Aj=j„-j a , &j=j T -j F 

k nodal diameter index 

L number of samples per blade passage 

/ counter relative to blade passage 

M moment (eq. (152)) 

m blade counter, increases opposite to direction of 
rotation 

m r blade counter, increases in direction of rotation 
N number of blades 

n denotes harmonics (eq. (43)) 

P unsteady pressure 

AP suction-side minus pressure-side differential 

pressure 

Q total number of revolutions 

q revolution counter 

R radius (fig. 5) 

r distance from torsional axis along blade chord 

(fig. 7) 
t time 

u frequency variable in E * s (eqs. (14) and (49)) 

V relative velocity 

W aerodynamic work (eq. (154)) 

x distance coordinate (fig. 5) 

hx deflection along x (fig. 5) 

y distance coordinate (fig. 5) 

z argument of the gif function 

a physical variable associated with a blade 

/3 defined preceding eq. (155) 

T defined by eq. (93), (103), (132), or (139) 

7 phase angle of G 

7 phase angle defined by eq. (129) 

e small quantity (eq. (151)) 

f distance from leading edge along blade chord 

7} distance coordinate (fig. 5) 

0 deflection angle (fig. 7) 

k nodal diameter index with range given by eq. (14) 
A defines limits on n (eq. (43)) 

X negative of lower limit on k (eqs. (63) and (48)) 

/x integer (eqs. (10) and (37)) 

v integer (eq. (36)) 

£ distance coordinate (fig. 5) 

p density 

a inter blade phase angle (eq. (1)) 


$ angular coordinate in rotating frame; increases 
in direction opposite to rotation 
angular coordinate in rotating frame; increases 
in direction of rotation 
<p phase angle of a (eq. (7)) 

X position of torsional axis as a fraction of chord 
from leading edge 
^ phase angle 

$ phase angle defined by eq. (41) 

go circular frequency 

aj r rotor angular velocity 

Subscripts: 

a relates to set a (eq. (17)) 

a subscripted variable is incorrectly associated with 
set a 
F false 

/ flutter 

/ imaginary part of complex variable 

k relates to k^ nodal diameter 

/ relates to /th point 

M relates to moment 

n relates to n th harmonic 

AP relates to differential pressure 

s steady state 

T true 

t total value (sum of steady and unsteady 
components) 

u relates to set u (eq. (17)) 

u subscripted variable is incorrectly associated with 

set u 

dx relates to bx 

0 relates to 0 

k relates to *th nodal diameter 

fi relates to p 

0 blade number zero 

Special symbols: 

S S is complex 

S only real part of complex variable S has a 
meaning 

S* complex conjugate of S 

S complex amplitude in real spectrum related to 

complex amplitude S in complex spectrum 
S' variable associated with second transducer; 

corresponds to S of first transducer 
AS where S~l,m,v, or $; denotes displacement 

increment of second transducer relative to first 
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Properties of Stationary Spectra 

Blade Sampling in a Stationary Reference Frame 

Since flutter is a self-excited vibration, it is generally 
associated with one of the natural modes of the bladed 
disk assembly. Usually therefore, the ratio of flutter 
frequency to rotational frequency (i.e., flutter frequency 
expressed in engine orders, or E 9 s) is nonintegral. It is 
because of this property that stationary sensors can be 
used to detect flutter. 

All blades in flutter vibrate at the same frequency with 
a well-defined phase relationship between blades. In the 
simplest case the phase angle between the neighboring 
blades (the interblade phase angle) is constant and the 
vibratory amplitude around the rotor is uniform. For an 
A-blade rotor the interblade phase angle a is restricted to 
N values given by 

a= 2 -^ k=0 N - 1 (1) 

This relationship follows from the periodicity condition 
that must be satisfied when completing a full circle 
around the rotor. The stability of the rotor system 
strongly depends on the value of a. This follows from 
both theoretical analyses and experimental observations. 
The index k also represents the number of diametral 
nodes (or nodal diameters) around the rotor or equiva- 
lently the number of complete cycles around the rotor at 
any given time. 

In terms of the wave equation, therefore, flutter can be 
described in a rotating frame of reference by 

a(k,$ n t) =A k cos(o)ft — $ r k + \f/ k ) (2) 

where a represents a physical quantity such as the blade 
displacement or the pressure on the blade surface; A k is 
the amplitude; a y is flutter frequency; is the angular 
coordinate, positive in the direction of rotation; and ^*is 
the phase angle corresponding to $ r = 0 and time r = 0. In 
this section sampling only at the blade locations will be 
considered, regardless of whether or not a physical 
variable is also defined in the space between the blades 
(i.e., static pressure). Thus, $ r can take on values only at 
locations corresponding to the blade positions, and 

* r =?£m r +B, m r =l, . . .,N (3) 

where m r is the blade number, increasing in the direction 
of rotation, and B is a constant. Therefore, 

a{m r ,t) =A k cos(u f t-^m r k-kB + ^^\ (4) 


In general, the blade amplitudes and interblade phase 
angles are not necessarily uniform around the rotor. This 
occurs, for example, when the single-blade natural fre- 
quencies are not all identical; that is, when the rotor 
blades are mistuned (ref. 6). Then a is given by a 
summation of possibly all AT waves; that is, 

N - 1 / 2tt \ 

a(m r ,t) = ^A k cos yuft- -^m r k-kB + ^ k J (5) 

These waves move with respect to the rotor. To find their 
speed, equation (4) is differenced (in the discrete domain 
this operation replaces the differentiation) with respect to 
time, holding a and therefore the argument of the cosine 
function fixed. That is, 

^(uft-^ntrk-kB + t^O 
or 

2?r A m r __ c cj 
~N At " k 

Thus the angular velocity of the k * h wave is inversely 
proportional to k and is positive in the direction of 
rotation for k<N/ 2. When k>N/2, the wave has the 
appearance of moving opposite to the direction of 
rotation. For this reason it is frequently more convenient 
to let -N/2<k<N/2~l for even N, and 
- (AT- l)/2 <A:< (N- l)/2 for odd N; then the sign of k 
indicates directly the sense of the wave propagation. For 
k> 0, the wave is forward, propagating in the direction of 
rotation; and for k < 0 the wave is backward, propagating 
opposite to rotor rotation. In the subsequent derivations 
different limits on k will be selected depending on a 
particular application. However, in general expressions 
the original limits on k, 0<k<A r - 1, will be retained 
since then no distinction between even N or odd N is 
necessary. 

In terms of complex exponential notation, equation (5) 
becomes 

a = a m «*•/' = A k eH* f t-2™ r k/N-kB) ( 6 ) 

where a denotes a vector only the real part of which has a 
physical meaning and 

|a J e^r, A k = \A k \ e^=A k (7) 

From the point of view of a stationary observer t and m r 
are not independent since 
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w r / M = -$ r + q2ir= - — m r —B + q2% 

q = 0 Q-l (8) 

where u r is the rotor angular velocity, q is the revolution 
counter, n is an integer, and t=t /1 is discrete. Let now 

B= -2ir, m=N-m n m=0, . . ., N- 1 

then 

2ir _ 2tt . . n 2 7r .a, 

«#■*,» = jj m +q2*= (m + qN) = jjix ( 9 ) 

where 

fi=m + qN /x = 0, . . QN- 1 (10) 

The blades are therefore renumbered so that for ?>0, 
m >0. Thus, in a stationary reference frame 

Oi = a m e'( w / w r) 2irn/N-_ Afc + ^ +A:m]2ir/^/ 

^ A: = 0 

TV-1 _ 

_ ^ ^/'(a;//w r + A:) (m + ^7V r )2x/7V 

A: = 0 

= ^ A* ^/ a r+k)2irfi/N (11) 

k = 0 

where in the second line use was made of the periodicity 
of the exponential function. Note also that m in a m could 
be expressed in terms of p by using the remainder 
function 

m=tx mod N=tx — N gif 

where gif(z) is defined as the greatest integer that is less 
than or equal to z and z is real. (Note that although in this 
equation p and N are integers, the mod function is de- 
fined in the same way when these two variables are real.) 

The real part of the left side of equation (11) represents 
a physical variable a that is being measured, and the real 
part of the right side can be viewed as a representation of 
the output of a digital spectral analyzer. In the last line of 
equation (11) the term in parentheses in the exponential 
can be interpreted as frequency in units of rotational 
frequency and 2tt/N as a sampling interval in radians. 
This representation was found to be convenient and 
therefore it will be retained in subsequent derivations. 

It follows from equation (11) that for each nodal 
diameter k there will be a separate frequency line in the 


spectrum of a stationary sensor. This frequency, ex- 
pressed in engine orders, is shifted away from flutter fre- 
quency by k. Because of this shift it is impossible to 
determine the true flutter frequency from the spectrum. 
However, as will be shown later, flutter frequency can be 
determined by using two stationary sensors mounted at 
the same axial location but at different circumferential 
locations. Also, in practice, the frequency can be deter- 
mined from blade-mounted transducers such as strain 
gages (which are frequently used with stationary sensors). 

Having interpreted equation (11) one can leave out the 
common terms on the two sides in the first line of the 
equation. This results in 

N-\ _ 

E A k ei2irkm/N (12) 

k = 0 

Except for the omission of the multiplying factor 1/N, 
this expression can be recognized (ref. 7) to be the inverse 
discrete Fourier transform (DFT-1). Therefore the A k s 
are given in terms of a DFT as follows 

A k =^ E a m e-nTkm/N (13) 

” m=0 

where the factor 1/Nis necessary because it was omitted 
in equation (12). 

Provided that the frequency shift w//« r in equation 
(11) is known, the amplitudes and phases for each mode 
can be easily determined. The details of this procedure 
include consideration of frequency aliasing and are 
presented in the next section along with the interpretation 
of individual blade spectra that result when m in equation 
(11) is held fixed during sampling. 

Interpretation of Folded Spectra 

Since m represents a spatial coordinate, equation (12) 
can be interpreted to give a spatial distribution of 
amplitude |a m | and the instantaneous temporal phase 
<P m . Because a m is complex, all N coefficients A k are 
independent. On the other hand, for N blades only N 
samples per revolution are available. Recalling that the 
number of samples per revolution is equal to the 
sampling frequency expressed in units of rotational 
frequency, from sampling theorem (see ref. 7, e.g.), it 
follows that the folding frequency is N/2. This implies 
that only N/2 coefficients for even TV and either (TV+ 1)/2 
(if ( W/ /av)mod l<l/2) or (N- l)/2 (if («/w r )mod 1 
>1/2) for odd N could be determined because the 
frequency associated with the remaining coefficients will 
be folded. It is shown later, however, that folded 
frequencies can be predicted so that except for the cases 
(wy/w r )mod 1 = 1/2, or 0, all N coefficients can be 
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obtained from the spectrum. As already mentioned, 
cy/oy is assumed to be known. 

To be able to determine coefficients in equation (11) 
from a spectrum, it is convenient to express this equation 
in terms of an integer k and a real variable u K9 which are 
defined by 


Note, however, that in equation (17), u K in the exponent 
has the same range as in equation (15). The frequency a K 
can be expressed in terms of u K (of eq. (15)) as 

a K =N-u K =N-(j£+?) ~<u k <N (20) 


0) r U) r 


► 


U.= ^ + K 


«r 


0<U k <NJ 


Equation (11) then becomes 


a(ji) = ^ A u, e* u ' 2rii/N 0 <u k <N 

K 


(14) 


(15) 


Note that complex amplitudes associated with folded fre- 
quencies, as shown in equation (17), must be conjugated. 
Complex amplitudes of equation (15) are given by 


A Uk ^(^u k 0 <u K < ^ 




^4 w — &N— w . 


N 


<u k <N 


( 21 ) 


They can be relabeled in terms of k as follows 


The u K is used in the subscript of the complex amplitude 
A instead of a more logical subscript k , because, to begin 
with, k’s are unknown and u ’ s are determined first from 
the spectrum. Since oy/oy is assumed to be known and k is 
an integer, it follows from equation (14) that 


u K mod 1 = ^mod 1 (16) 

0) r 

This could be used to identify peaks associated with the 
set u K without the knowledge of k’s. However, because, as 
already mentioned, the data are undersampled, some fre- 
quencies associated with ay/ ay will be folded. In general 
therefore, the stationary spectrum of a measured variable 
a can be represented as 


^4 k ^4 Uk _ 


( 22 ) 


and in terms of k (eq. (12)) as 
A[ k) =A^ N k<0 


(23) 


A^—A^ x>0 

where the ranges of positive and negative k’s, assuming 
that a )f/u r <N, are 

0<k<N-^ - ^<k<0 (24) 

0) r Oy 


d(#i) = £ ( a u K e' u * 2 *^ /N + a* K e iu ^ /N ) 


By substituting the relation between k and k implied by 
equation (23), it follows that 


where G denotes a complex amplitude and a K denotes a 
folded frequency set distinguishable by the property 

a K mod 1 = 1- (cy/oy) mod 1 (18) 

Because there can be at most N complex amplitudes G 
for a given k , only one (either GL Uk or G a ) can be 
present in equation (17). The ranges for frequencies 
associated with G’s are 


( 19 ) 


„ N 

o<u K <- y> 


c Of N cx)f 

£<K< — J- 


<K< . 

ay 2 


2 2 ay ay 


0<k<N or 0 1 

as assumed in equation (12). 

Note that because of equations (16) and (18), when 
(oy/ay)mod 1 = 1/2, the two sets of frequencies u K and a K 
cannot be separated, and consequently, without addi- 
tional input the k’s (and associated A k’s) will remain un- 
determined. However, this situation, just like the case 
where cy/ay happens to be an integer, can be avoided in 
practice by slightly shifting ay from the particular value. 

Once all A^’s are determined from a stationary spec- 
trum, the inverse DFT (eq. (12)), can be used to obtain 

Note that the limits for the two ranges of k given by 
equation (19) are not inclusive. Although these limits 
sufficed for the preceding derivation, for later reference 
it is useful to derive the precise limits on k’s. 
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Using the mod function in equations (19) gives 

_ (^1 — i^mod A (-y- - 

\ u r u r ) \ 2 «/■/ 

-(^_2f modI ) (25, 

_ /w£_f^ mod A (26) 

\w r C3 r / 

For even iY, therefore 

_ ( a- A SKuS " - 1 _ (a - a mod A 

\«r *>r / 2 \<*V o r / 

(27) 

(^-^mod l) </c a <(AT-l) 

2 \ w r w r / 

— — ^mod A (28) 

V<*V w r / 

For odd TV the distinction must be made whether 
(wy/co,) mod 1 is less than or greater than 1/2. Thus for 
(coy/oj r )mod 1 < 1/2 

- - ^mod A </c u <^-^ 

W «r / 2 

_^_^ mod A (29) 

V^r “r / 

- ^mod A <K a <(N- 1) 

2 \u r u r ) 

-(^Z — ^mod A (30) 
\«V 03 r / 

and for (wy/co^mod 1 > 1/2 

- — ^mod A ~ — ^mod A 

\<*V co r / 2 \co r 0 ) r ) 

(31) 

-(ZZ - ^mod A <K a < (AT- 1) 

2 V“V “V / 

-(^_“r mod A (32) 

\«V w r / 


Consider now sampling only one particular blade; that 
is, m during sampling is considered fixed. (Actually, 
samples can be collected from all blades as they pass the 
sensor. However, the samples should later be regrouped 
into AT groups so that any one group will contain samples 
from one particular blade only.) Then, 

&(m,q) =A m e' ( “/ / “r)2T 9 g=0> . . q_ ] ( 33 ) 

To avoid aliasing, coy*/ co r must be less than 1/2. For a 
practical range of c oy* and co r this case can be ruled out, 
and therefore a single blade spectrum will always contain 
aliased flutter frequencies. 

If w= (coy/oo^mod 1 < 1/2, the flutter frequency will be 
folded an even number of times and hence 

A m =(a u ) m ( 34 ) 

If « = (coy*/co r )mod 1 > 1/2 in the spectrum, the aliased 
frequency a , a= 1 -(coy/co r )mod 1, 0<a<l/2, will be 
folded an odd number of times and 

A m =(Q a )‘ m (35) 

Note that the overall (i.e., N samples per revolution) 
stationary spectrum is folded only once about N/2 
engine-order frequency, whereas a single blade spectrum 
is usually folded several times about 1/2 engine-order 
frequency. This implies, first, poorer peak separation in 
the single blade spectrum than in the overall spectrum. 
Second, because most of the noise in measurement 
systems falls in the low-frequency range (in this particular 
case of the order of 1/2 E and below), the signal-to-noise 
ratio of a single blade spectrum will frequently be below 
that of the overall spectrum. And third, the overall 
spectrum essentially gives a modal (i.e., in terms of 
spatial modes associated with each k) expansion of a 
physical variable, which is in general of greater interest 
than the individual blade characteristics described by a 
single blade spectrum. For these reasons most quantita- 
tive analyses will be presented in terms of the overall 
spectra. (An exception to this would arise, for example, 
during flutter monitoring, when one blade exhibits par- 
ticularly large deflections. Its displacement spectrum is 
then desirable, and it is most expediently obtained by 
sampling only the particular blade.) 

Figures 1 and 2 (ref. 4) illustrate the displacement and 
differential pressure spectra obtained by sampling N 
(N= 38) times per revolution. The differential pressure 
represents the difference between the suction and the 
pressure sides of the blade and is nondimensionalized by 
using the inlet dynamic head. The peaks in these spectra 
are identified by nodal diameters k, where the range of k 
can be obtained from equation (14) by substituting TV and 
the true flutter frequency coy/ co r = 8.45; that is, 
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1-0 r 


.8 h 


= .6 


E 

a> 

fO ^ ij. 


-5 


k = 0 <8. M*5) 


* = 20 
(9.55) 




0 2 


J 

8 10 12 1H- 16 18 20 


Frequency, engine orders 

Figure 1 Displacement amplitude spectrum; port A; 1 unit of displacement = 3.267 X 10~ 3 radian. Average of four FFT’s. 



Figure 2. -Differential pressure amplitude spectrum; port 1; amplitude nondimensionalized by using inlet dynamic head. Average of four FFT’s. 


-8<k<29 

Although other less significant nodal diameters associ- 
ated with (oy/co r )mod 1=0.45 were not marked in the 
figures, they can be easily identified by simply counting 
peaks to the left and right of the zero peak. The most no- 
ticeable response associated with 1 — (a>//ci> r )mod 1 =0.55 
is the peak at ay/co r = 9.55 Oust to the right of the #c=l 
peak at ay/o> r = 9.45). From equation (20) it follows that 
the corresponding k for this peak is 20. 

One of the essential steps in obtaining a pressure- 
difference spectrum, such as in figure 2, is determining 
the blade pressure and suction surface locations. In other 
words, one needs to establish a criterion by means of 
which an JV-ordered pair of data points representing 


blade surface pressures can be selected from a continuous 
stream of data points that encompasses the pressure field 
for the whole circumference of the rotor. This was 
accomplished by using the steep pressure gradient region 
in the steady-state distribution (ref. 3). 

Complete Pressure-Field Sampling 

In this section a theoretical expression for the complete 
pressure spectrum will be derived in a form that is most 
closely related to the measured spectrum obtained by 
using a fast Fourier transform (FFT) analyzer. Once the 
interpretation of the measured spectrum (discussed in the 
next section) is completed, the theoretical expression can 
always be recast in a more compact form, as shown later. 
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Since pressure is continuous, the number of samples 
per revolution is not necessarily limited to the number of 
blades. Therefore, assume that there are L sampling 
intervals, /=0, . . ., L-l per blade passage m, where 
m = 0, . . N-l. The spatial dependence of the 
unsteady pressure can now be described in terms of 

v = l+mL y = 0,. . .,LN-1 (36) 

However, to describe the time dependence, another vari- 
able fi needs to be introduced. 

H = v + qLN ^ = 0, . . .,QLN-l (37) 

Note that for L = l, /x is consistent with the definition 
given by equation (10). 

Assume that for one particular blade passage (m = 0), 
and for a given nodal diameter k, the pressure amplitude 
distribution is given in terms of the discrete set / by f k j. 
Similarly, let the discrete temporal phase distribution be 
given by \}/ k j. Then in the rotating frame of reference the 
pressure field can be expressed as 

P(k,l,t) =fk,i cos(o> f t+$ k j) 

and the instantaneous pressure distribution as 

P(k,l t 0) =f k j cos 4r kJ 

For an arbitrary blade passage m, P becomes 

P(k,l,m,t) =/*,/ cos (a) f t + + yl'kjj (38) 

(Recall that m increases opposite to the direction of 
rotation.) In a stationary reference system, in terms of a 
complex exponential notation, P becomes 

P(k,l,m,n)= \f kJ \ e ih»f/o>rn*>*/LN+ 2 irkm/N+t k ' l ] (39) 

where 



Adding and subtracting lirkl/LN, adding 2t kq in the 
exponent, and then using equations (36) and (37) result in 

P(k J i,fi)~ \f k /| + k) 2xfifLN-2rkl/LN + t kJ ] (40) 

Letting 

*kj=hj+™ (41) 

in equation (40) gives 


P(k,l,n) = \f kJ \ efa eH»//ur+k)2w/LN ( 42 ) 

The first term in this equation can be expanded as 
follows: 

fk,l= \fk,l\ ih ‘‘= L E A Fk.n e^ln/L (43) 

n= -A 

where A=L/2 for even L and A= (L~ l)/2 for odd L 
and F k n 9 s are given in terms of a complex DFT by 

1 L-\ __ 

F k.n=~r ^ fk,ie~ nrln/L ( 44 ) 

If n is restricted to be positive, equation (43) becomes 
- n*(L- 1)/2 

fk,l= £ ( F k,-n e-n*i*/L+F k , n ei^ln/L) 

n= 1 

+ Pk,0 + Pk,-L/2 (45) 

where the upper limit on n is (L - 1)/2 for odd L and 

L/2 - 1 for even L , and the last term is included only 
when L is even. Because the inverse DFT is periodic 
beyond the sampling interval, /x can be substituted for / in 
this equation. Performing this substitution and then sub- 
stituting the resulting expression in equation (42) gives 

n*(L- 1)/2 

P(k,fi)= £ (p k _ n e i ( a f /u r+k-nN)2* f L/LN 

n- 1 

+F k n & ( w // w r + £ + nN) +F k q e l ( ^f/^r + k ) Ir^ilLN 

4. F k _ L/2 j{<*f/(j3 r +k-LN/2)2Tn/LN (45) 

If, in addition, the summation over all possible k is per- 
formed, equation (46) becomes 

N- 1 Vn< t {L-\)/2 / 

P( M )= E E ( F k,-n e‘ ( “f / *r +k - nt *)2r»/LN 

k = 0 [_ n = 1 

+F k „ g‘( a // u r+k+nN)2rn/Lf^ +F k Q gi(<>>s/Wr+k)2x)i/LN 

+ F k _ L/2 e*( w f /a r +k-LN/2)2Tii/LN\ ( 47 ) 


where 2'k^i/LNxs proportional to time. In the next section 
this equation is used to obtain k’s and n’s from the mea- 
sured folded spectrum. 
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Interpretation of a Complete Pressure Spectrum 

To interpret a complete pressure spectrum, it is again 
convenient to introduce k instead of k, where as before 
(eq. (14)) 


, 1 <k<N 

co r co r 


or 


(48) 


- ^mod \ ) </c<7V— 1 - — ^rnod l\ 

\0) r 0) r / \ U/- / 

Denoting frequency terms by 

«K,±/t=^ + K±nN (49) 

OJ r 

equation (47) can be expressed as 


The main purpose of this equation is to define the 
notation and to illustrate the procedure that follows. 
That is, it is attempted to^ associate the aliased ^F’s with 
u Kt -Vs and the unaliased £Fs with u KtfJ ’ s. As proved later, 
this is correct except for n= (L- 1)/2, and therefore this 
term was singled out in equation (51). Note that the 
ranges for u K n and a Kifl derived below refer to the sub- 
scripts of Ts and that u Kt±n ’ s in the exponents are the 
same in equations (50) and (51). 

First consider the set of frequencies u K>n . Substituting 
the limits on k from equation (48) into equation (49) gives 


nN<u Kn < (n + l)N n = 0, . . or 

or 


‘ 

> (52) 


nN+ ^rnod 1 <«*„< (n + l)N - 1 + ^mod 1 

O), w. 


K L n = 1 

4. F Ju^^ii/Lbh . F Ji4 K ' 0 2irii/LN 

T r t*K,n K c 

+ F“'.-U2 e ‘ UK - ~ i/22,r ' i/LA 'J (50) 

The u Kn ’ s in the exponent and in the subscripts of F’ s 
are used because, to start with, k ' s and n’s are unknown 
and these frequencies are the most readily obtained 
quantities from the measured spectrum. However, once 
k’s and n’s are determined, it is more convenient to 
denote the complex amplitude as F Kifl . 

Because the sampling rate is LN per revolution, folding 
is expected in the measured spectrum for u Kttl >LN/ 2. 
Therefore for this spectrum equation (50) is replaced by 

[ n<{L-X)/2 , 
n- 1 


In addition, the subscript u K>n in equation (51) (but not 
that in (50)) is restricted by definition to 


For odd L and n=(JL- l)/2, it is this condition that 
limits the subscript u K>n at the upper limit in equation 
(52). Except for this case, it follows that, by subdividing 
the measured spectrum into segments given by the in- 
equalities (52), each segment will correspond to con- 
secutive values of n, n = 0, . . ., L/2 - 1 or (L - 3)/2, and 
all k frequencies associated with each n will be contained 
entirely within each segment n in the order described by 
inequalities (48). Explicitly, for each n , k’s associated 
with each u Kt „ can be determined from 

K = u Kt „-~£-nN (53) 

where 

u K n mod 1 = ^mod 1 (54) 

Cx) r 


_|_ cy e iu Kin 2TCix/LN\ + g: ^u^irfi/LN 

U K,n / 0 


+ *Kui e< - i/22TM/ V + E $«, 


J K < N/2 - U3f/(j} r 


%{L~ l)/2 


is used to identify the peaks that are associated with the u 
set in the measured spectrum. For odd JL, n= (L- 1)/2, 
and u KtU <LN/2 , upon substitution in equation (53), 



a )f 






+ 


£ 

K>N/2~<j}f/b) r 


C V* 

° t*K.(L -l)/2 


e * u K, [L - 1 )/2 2 Xfi/LN 


(51) 


For these k’s and n’ s, the F Uk ’s (or alternatively F K>n ’s) 
are given by 

Fu K . n =FK,n = 5u Kin 
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For the particular case n= (L- 1)/2 and K>N/2-o)f/<jo n When L is odd, it follows from equation (60) that for 

it is shown below that all n’s 


« = (-«*,«) - 7^ + 
co r 


L + 1 L-l KT LN ... , 

~t~ n —*- N<a K,n<-r < 55 > 


a n,n< ~~ 2 ~^ 


where 


Therefore, the range 


a Kfn mod 1 = 1- °^mod 1 
c c r 


(56) 


L-l 

2 


M<a Ki „< 


LN 

2 


and 



Consider now u Kt _ n frequencies for which 


which is bound by the aliasing condition at the upper end, 
must be associated with u Kytl for n=(L- 1)/2 and 
K>N/2 — Wf/(jo r . (Recall that u Kyn ' s were previously 
defined only for k<N/2 — ay/av for this n.) Therefore 
substituting 


- nN< u K> - n < - (n - 1)TV 



or 


L- 1 
2 


or V (58) 



Since all of these frequencies are negative, they will be 
associated in the measured spectrum with the set a; that 
is, 

u k, — n~~ a K,n (^9) 

and 


u K,n—LN—a K>n n= (L- 1)/2 
in equation (53) gives 


K=LN-a lin 



L-l 


N 


= (~a K , n ) - 



L + 1 
2 


N 


which is equation (55). Recognizing the fact that a Kyn f s 
are folded frequencies, equations (56) and (57) follow 
directly. 

The particularly interesting case results when L = 1 . 
Then equation (55) reduces to (20). Thus the results of the 
section Interpretation of Folded Spectra, obtained for a 
sampling rate of TV per revolution, are recovered. 

For practical reasons it is also desirable to express the 
results in terms of F^s where A: = 0, . . ., TV- 1 and 
n^- 0, . . ., TV- 1. The necessary transformation is given 
for TV— A<A:<TV- 1 and -X</c< - 1 by 


(n-l)N<a K '„< nN n=l or 

or (60) 

(n - 1)TV+ 1 - ^mod 1 <a K n <«TV— ^mod 1 
ay K,n ay 

These frequencies can be recognized in the spectrum from 
the property (56). For each value of n. 


K = (-a K , n )-jf+nN 

w r 

and 

F u = =F K __ n = $* a 

u k, -/I n u K n 


(61) 


(62) 




0<«^<L-2 — A 


1</2*<L-1-A 


Fi k X 


= f { k ) — /r(*) 

n k * k— N,n k — L + 1 **—n K 


fL-l-A 
1 — — A 


Fi k X- l=^M 0 = F $ (63a) 

and for 0 <£<TV— X- 1 and 0</c<TV- 1 - X by 
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p{k) — pM 

1 k,n k 1 k,n k “~ K,n K 


0<n k <L - 1 - A 


0<n<L - 1 - A 


p(k) —pM p{K) 

k,n k 1 k,n k -L * K,-n K 


L — A<n k <L- 1 


1 A 


(63b) 


where A= (L- l)/2 for odd L and A—L/2 for even L 
and - X denotes the lower limit on k in equation (48). The 
k and k in the superscripts on the F 9 s and in the subscripts 
on the n 9 s were used to clearly distinguish between the 
two families of F 9 s and «’s. Note that for L = 1 equations 
(63) reduce to 


which again involves the use ofthe inverse DFT. From 
the real and imaginary parts of P the pressure amplitude 
and the instantaneous temporal phase distributions in a 
rotating frame of reference can then be constructed. 

Equation (64) will now be used as a starting point to 
prove equations (63). To simplify the derivation, it will be 
assumed at first that n k and n K have the same range, that 
is, 


n k — Q, . . N— 1 n K = 0, . . ., N- 1 

Substituting equation (43) with A = 0 into equation (64) 
gives 


0 ) = 


N- l 



L ~ 1 __ 

^ Fjf^ etenlnk/L e i2irkl/LN ei2ickm/N 
n k = 0 


= N-\*KSN-l —\<k< — 1 

0<At<jV-X-1 0<k<N-\- 1 Since 


(65) 


Equations (63) will be derived later. 

Thus, all n and k can be determined from the pressure 
spectrum, again, provided that (oy/co r ) mod 1 is not 1/2 
or 0. Note that n for both families, u and a, increases 
from left to right in the spectrum. However, for aliased 
frequencies within each n, k varies from right to left. 
Since in the spectrum in each n segment there are two 
different values of n, one associated with set u and the 
other with set a , to distinguish between the two, one 
could associate the negative n 9 s with set a. The only 
exception to this rule arises when L is odd, and 
k > N/2 — co r . The aliased frequency in this range is 

associated with n— (L- l)/2. (Negative n 9 s could also be 
retained in the expressions associated with set a , provided 
that absolute signs are included with each n.) 

Once the F k n 9 s are determined, as shown by equation 
(43), they can be used as inputs to the complex DFT-i 
algorithm to obtain l for each k. From real and 
imaginary parts of f k ( the distribution of pressure 
amplitude \fk,i\ and the instantaneous temporal phase 
\p k j (eq. (41)) can then be obtained for each k. This phase 
distribution really applies only for blade passage m = 0. 
For any other m the entire phase angle in rotating frame 
of reference (eq. (38)) is 


k = K + N — A < — 1 N-\<k<N- 1 

k = K 0 </c 1 — X 0<&<7V— 1 — X 

only the following part of the summation with respect to 
k from equation (65) need be examined: 


N- 1 L- 1 

Y y pi k ) e i2ivln k /L e i2irki/LN eilirkm/N 

k=N—\ nft 0 k ' nk 

where -X denotes the lower limit on k in equation (48). 
After substituting k = K + N and using the periodicity 
property of the exponential function, the preceding 
expression becomes 



y p(k) J2wl(n k + \)/L 
r K+N,n k e 

n k = 0 


gi2i tkI/LN eil-KKni/N 


Now substituting n K =n k + 1 gives 



rt K =l 


Fi%,n K - 1 


gilTrlnjL ^2 tt kI/LN (> ilirKm/N 


'l / k,i+ 


2-irkm 

sr 


= $k,l + 


2irk(l+mL) 

LN 



i-i 

2 ^ /?(*) e >1nlnJL e i2 VK l/LN e i2mm/N 
n K = 0 K, " k 


( 66 ) 


The instantaneous pressure distribution around the 
rotor can be obtained as follows 

P(l,m,0)= N f} \f kJ \ eilhl+l*k(l+mL)/LN] (64) 

k = 0 


The right side of equation (66) corresponds to the part of 
the summation for negative k that would be obtained if k 
were used instead of k in equation (65). Therefore term- 
by-term comparison of the two sides of equation (66) 
gives 
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p(K) w ?(k) = /?(£) 

x n,n K * K+N,n K _ i J k,n k 


1 <n K <L- 1 


►-X<*< - 1 


In the remaining range of k 


pM = /r(fc) 
* K,n K A k,n K 


0<n <L — 1 X — 1 


These equations can also be expressed in terms of k and 
rife as follows 


^-1=^-V,0 


0<n k 



N-\</c<N - 1 



0<« t <L-l 0<A:<N-\-l 


(67) 


It is now necessary to derive the transformation 
between the Fj*% and F^J^, where the range of A: must be 
the same for both F’s and 


0 <n <L - 1 ~ A<n^<L - 1 - A 


In the derivation, the range of k is assumed to be X: = 0, 
. . N- 1; however, the results do not depend on this 
choice and any other range could have been selected. 

When we use the preceding notation, equation (65) 
becomes 


_ N — 1 L ~1 

P(l,m,Q)= ^ 2 F^ eil'Kln fL gUj-kl/LN eil-wkm /N 

k = 0 /i = 0 


The part of the summation for L — A <n <L - 1 is 


The right side of this equation was obtained from the 
expression for P(l,m, 0), written in terms of by 

including only the part of the summation corresponding 
to negative /* A . Term-by-term comparison of the two 
sides of the preceding equation then gives 




-A</z a <-1 


or 


F iVn= F k,n-L= F &l £-A<«<L-l > (68) 

0<«<L — A- 1 


where the last equation was added to complete the range 
on n. 

Equations (63) result by consecutive application of 
equations (67) and (68) with an additional slight modi- 
fication resulting from the restriction that n be positive. 
Thus the negative sign must be carried into the second 
subscript of F. 

Figure 3 illustrates the complete unsteady pressure 
spectrum for the same flutter point as for the spectra in 
figures 1 and 2. Hence ay/co r =8.45 and Af=38. For the 
set u frequencies 


%nod 1=0.45 
0) r 

which identifies the associated peaks in the spectrum. The 
spectrum can now be subdivided into segments each 38F 
long. For the frequency set u, the first segment corre- 
sponds to n = 0 and the last (which is only partly in- 
cluded) ton = 6, Equation (53) is now used to determine k 
associated with each u frequency marked in the particular 
n segment. The resulting k s and n ’ s were also marked in 
the figure. 

For the set a 




eilrln/L gilrkl/LN eilickm/N 


Substituting n=n\+L in the preceding expression and 
then using periodicity of the exponential function result 
in 


N - 1 -1 

E E 

k=0 „ A =- A 


F i n n K +L e i2vl " A/L e^xkl/LN e i2xkm/N 


N-l 




p(n A ) JhlnJL e i2ikl/LN e ilxkmlN 
K > n a 


l-^mod 1=0.55 

u r 

and the 38 E segments in the figure were numbered 
beginning with n = 1 and ending with n = 7 (which is 
incomplete). Equation (61) was then used to determine 
k’s associated with each a frequency marked in the 
particular n segment. Once k’s were determined for this 
set, they were negated to distinguish them from positive 
n’s associated with the set u. The resulting k’s and n’s 
were then marked in the figure. (Note that one can 
introduce negative n’s with this set of frequencies from 
the start provided that n is replaced with \n \ in eq. (61).) 

Another example of a complete pressure-field 
spectrum obtained for an earlier build of this fan rotor 
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Frequency, engine orders 


n = 0 
FRQ, <E> 

K 

r->= - 1 
FRQ, (E) 

K 

T-, = 1 

FRQ. < E ) 

K 

r*- -2 
FRQ. CE) 

^ = 2 
« FRQ. (E) 

K 

4-. 4-5 

-4- 

25 o 55 

4- 

4-1 . 4-5 

-5 

6 4-o 55 

3 8 1.4-5 

-3 

5 o 4-5 

-3 

26 o 55 

3 

4-2 o 4-5 

-4- 

67.55 

0 87.4-5 

3 

6.4-5 

-2 

2 9 o 5 5 

0 

4-3 0 4-5 

-3 

70.55 

-3 


8 o 4-5 
1 1 o 4-5 

0 

3 

32 o 5 5 

-3 

4-9 . 4-5 

3 





(a) To 250 engine orders. 

(b) Zero to 100-engine-order portion of spectrum in fig. 3(a). 

Figure 3. - Power spectrum of leading-edge pressure transducer signal; engine A. 


(engine B) is given in figure 4. For this figure, 
o)//co r =7.19 and as before N=3S . The response at 
aj/u> r = 205.75 is not associated with flutter. 

Note that in figures 3 and 4 only the low-order positive 
and negative k’s and n 9 s are present; this justifies the 
particular choice of limits on k and n . In terms of these 
k 9 s and /i’s, equation (47) reduces to 


AT-l-X Z,-l-A_ 

= ^ ^ F kH e i( Ct3 f/° > r + K+ nN) 2 Tfi/LN 

k= — \ n= — A 

where - X denotes the lower limit on k given by equation 
(48) and A is defined as in equation (43). 

In these examples the total number of samples was 
QLN= 8192. Conveniently this number is a power of 2 so 



13 




Frequency, engine orders 
(a) To 250 engine orders. 

(b) Zero to 100-engine-order portion of spectrum in fig. 4(a). 

Figure 4. - Power spectrum of leading-edge pressure transducer signal; engine B. 

that the standard DFT algorithm could be used. Because This will also be the most frequent case in practice, 

the number of revolutions Q was 16, the sampling rate It likewise does not seem to be important that Q be an 

LTV was 512. Note, however, that for 7V = 38 this results in integer. If so desired, however, both Q and L can always 

a nonintegral value of L (or more precisely an average L). be chosen to be integers by using Singleton’s algorithm 

Strictly speaking the preceding formulas apply only when (ref. 8). 

L is an integer. However, note that L enters the formulas Another practical consideration is the fact that in 

through the product term LN, which also determines the general ay/av will not coincide with any frequency line 

folding sampling rate, and either the L/2 or (Z,-l)/2 (i.e., will not be a bin-centered case). To evaluate the 

term in the summation over index n. In the examples, the complex coefficients from the spectrum, it is then best to 

sampling rate selected was sufficiently high so that the use the interpolation formulas that were originally 

terms associated with the tail end of the spectrum were derived by Donald C. Braun of Lewis. For completeness 

negligible and hence no difficulties were encountered. these formulas are included in the appendix. 
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Determination of Flutter Frequency 

The true flutter frequency can be determined from a 
cross-power spectrum between two pressure or two dis- 
placement transducers that are located at different cir- 
cumferential positions. The transducers must, however, 
correspond to the same axial position. 

To examine the spectral relationship between two 
transducers, an expression analogous to equation (46) or 
(50) but for an arbitrary spatial reference point is needed. 
By using this expression, an equation can be derived for 
the complex cross-power amplitude. From the cross- 
power phase angle and the known angular separation 
between the two transducers then, the particular nodal 
diameter (provided that the sampling rate is at least AO 
and flutter frequency can be found. Only one complex 
amplitude is needed from the measured cross-power 
spectrum. Therefore it is best to choose the maximum 
amplitude. As will be seen, the angular separation 
between the transducers can be selected to optimize the 
determination of frequency. 

Consider first the complete pressure- field spectrum. 
The original spatial reference in the rotating frame, m = 0 
and / = 0 at / = 0, will be retained for the first transducer. 
For the second transducer it will be assumed that, at t = 0, 
m = Am and /= A/ in the rotating frame of reference. For 
this transducer it is convenient therefore to introduce the 
following set of variables: 


in equation (71) results in 

P(k,l f ,v' ,t) ^ettuft + lrkv'/LN) e i2TrkAv/LN 


X E F k,rt ei^nAl/L e i2xnl'/L (72) 

n = — A 

Equation (72) can now be transformed into a stationary 
reference frame by using 

so that 

P(kJ' = +to'] 2 */LN e i2ickAv/LN 

L — \ — A __ 

X Ffc n e&irnAI Qi2TnV/L 

n— —A 

_ e i ( 0 ) f /ta r + k ) 2-Kfj. '/LN e j2ickAv/LN 
L- 1-A _ 

x E F k,n e' 2 ™ A/ e^nl’/L (73) 

n— —A 


m ' = m— Am 


/'-/-A/ 


v' — v~Av 


m' = 0, . . N—l; ^ 
0<A m<N— 1 

/'= 0, . . .,L- 1; 
0<A/<L-1 

v' =0, . . ., LN- 1; 
0<Av<UV-l 


(69) 


which vanish at f = 0, just like the corresponding un- 
primed variables for the first transducer. 

Equation (38), written for the rotating frame of refer- 
ence, can be expressed as 


P(k,l,m,t) = \f k> ,\ e'+ki e Ho> / i+2irkm/N) ( 70 ) 


(Note that both and equal zero when t is zero.) 
Since 

eilrnl'lL — qIIth (/' + m'L + q'NL ) /L = ' /L 

equation (73) can be expressed as 
Z, — 1 — A _ 

P(k,n') = e‘^^ v /LN £ F kn e i2*nAI/L 
n = -A 

X e i {o3//(j) r + k + nN ) 7 LN (14) 

Summing now over k results in 


by using the complex exponential notation. Introducing 
equations (41), (43), and (36) into this equation results in 

L-\-A __ 

p(k,l,v,t) = e ‘(<*f‘ + l*kv/LN) £ p k n e i2Hn/L ( 7 1) 

n— -A 


N-\ L- 1-A_ 

= T) e*2x*A v/LN V f k n eilrnAIfL 
k = 0 n~ —A ' 

^ pi { ay/a ) r + k + nN) Ixp 7 LN 


(75) 


Substituting now 
/=(/'+ A/) mod L 
v = (v ' + Av) mod(ZJV) 


Note that the expression for the first transducer can be 
recovered from equation (75) by setting A/=0, and 
Ay = 0. 

If n is restricted to positive values and if k (eq. (48)) and 
u K n (eq. (49)) are introduced in equation (75), it becomes 
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P(fi')= ^ e i2*K&v/LN 

K 



g — ilnnAl/L 


the corresponding amplitudes in the complex spectrum. 
Thus the cross power from the real spectrum for the set a 
is defined by 


x e iu *> - n2lCfi ' /LN + F„ 

U K,n 


e i2TnAl/L giu^Zru’/LNj 


&,,, -(®0 J KJ" 


+p Ju^irn'/LN+F e-'>A/ e‘<-L/22*nVLN I 

u k,0 C U K,—L/2 1 

[ n<(L- l)/2 / 

£ {K- 

n = 1 


ju Kt _ n 2’K\i'tLN 

+ F e iuK ’ n2ir t i /LN ') + F e iu K '02Tn'/LN 

U K,n J u k, 0 

+ —l/i e< ■ i/22 ™ /LN J (76) 


-<*0* 

= If I 2 e - '2^ ( KjAv/LN— n a Al/L ) 

\ u tc,-n I 

=o; > .,-|o^_.| = |S.„| A- (80) 

where again use was made of the relationship between 
F' and F implied by equation (76). Therefore 

( Av A/ \ 

Ka ZN~ n °2 ~) 

and 


In this equation primes denote the complex amplitudes 
for the second transducer. 

The complex cross-power amplitude between the first 
and second transducer can therefore be expressed as 

G K , ±n =\G K>±n \e i ^=(F' K>±n ) ( FU*)‘ 

= (K ±n enr^/LN e ±i2*nAI,L) (p K> ± „) * 

= \F k> ±„ I 2 e i^(KAv/LN±nAI/L) (77) 



where the integer^ n depends on k and n for a selected A / 
and Av. For set u, F Uk n = $ n and from equation (78) it 
follows, except for n=(L- 1)/2, odd L, that 






(83) 


where use was made of the relationship between F' and 
F implied by equation (76). Hence 


7 k, ±rt + 27T/ 


<,±n = 2ir^ 


kA^ nA/ \ 
LN ± ~L ) 


0— 7k, ±n<2ir 


and 


k = 


2JV 

Av 


( 7k, An 

V 2ir 



„A/ 


(78) 


(79) 


In principle, therefore, one can select some large peak in 
the complex cross-power spectrum, obtain y K> ± n , solve 
for k from equation (79), and then determine w//w r from 
equation (49) by using known k. In reality, however, 
because there are two sets of frequencies, u and a, in the 
measured spectrum, first a criterion must be developed 
that will permit one to differentiate between the two sets 
of frequencies. This is necessary because, as shown be- 
low, a different set of equations is associated with set a. 

Script letters will again be used to denote the complex 
amplitudes in the real spectrum to distinguish them from 


Since it is not known in advance whether the selected 
amplitude is associated with the set u or the set a, it is 
possible to associate it with the set u when in fact it is 
associated with the set a. For this particular case, 
subscript u will be used, with the dot indicating that the 
complex amplitude was incorrectly assigned to the set u. 
Thus, integer will be determined from the spectrum 
according to the rules formulated for the set u (in the 
section Interpretation of a Complete Pressure Spectrum), 
and subsequently k will be determined from equation 
(83). Hence 


— 


LN 

Av 




(84) 

i 


where integer n depends on k and n „ for a given A / and 
Av. In turn, c of/u r depends on k „ (eq. (53)); that is, 



Thus, in general (whether for frequency set u, u , or a) one 
must first guess the value of j to determine k and oy/co r . 
The resulting k must be an integer and the resulting ay/w r 
must be in the range 


16 


0<%-<N 

«r 


( 86 ) 


Substituting of from equation (89) in equation (91) 
then yields 


Thus j and k will have to be determined iteratively. 
However, in the case of an incorrect frequency assign- 
ment to the set u (i.e. , eqs. (84) and (85)), it may not be 
possible to determine physically plausible values for k and 

U)flO) n 

Because by assumption the set a is the correct one, 
there must be at least one set of values for K a , J Qk ^ and 
{<uf/G) r ) a such that the first two are integers and the third 
one satisfies equation (86). Since y tixn and y aKn are 
associated with the same complex amplitude 



/"A V A iK,n + M/L 

\u r ) a Am + Al/L 


Let now 


(92) 


p— + ^ 

A m+Al/L 


(93) 


With this definition, applying the mod function to the 
left and right sides of equation (92) gives 


Using then equations (82) and (84) and n a —n ^+ 1 in the 
above it follows that 


Av , riuAl 


Av 


A! . 


LN Ku+ ~ TN Kq+ ( - rlu+ L ja ^ 

and 

Al ^ LN ( . . x 

K,i ~ N Ap Ka+ ~Kv^ J <n Ja *J 

Denoting 

A/«c,n = Ju K „ ~Ja,, n 

equation (87) can be expressed as 


(87) 


( 88 ) 


(^)r d, ‘[["~(^)J* r ) modl 

#-(£>] mod 1 - T mod 1 1 


Then, if 


xmod 1 


•[[''(^)”° di ]' rmod, j 


xmod 1 


-K a +N 


^Un + Al/L 
Am+A//L 


(89) 


Note that it is only necessary to consider Am>0 and 
A/>0 since this can always be achieved by a particular 
choice of the transducer labeling. (This will also be 
assumed in all subsequent equations.) 

On the other hand, because u Kfn =a Ktn > using equation 
(85) and the following analogous equation for the set a 
(eq. (61)) 

(^0 =<-«*.«> -«c+n a N (90) 

results in 


T mod 1 < 1 


mod 1 


-(£>: 

mod 1 = 1- (jf'j mod 1 -T mod 1 


(94) 


and if 
T mod 1 > 1 


-09: 


mod 1 


mod 1 = 1- |r mod 1 — j^l — mod 1 J | 



+ K, i + n ti N= 

U 


(£)-, a + no N 


and 


=2 - mod 1 - r mod 1 (95) 

If T mod 1=0, referring to equations (89) and (93) it 
follows that /c„ will be an integer; that is, 

K t j mod 1=0 
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However, then from equation (94) 


and from equation (78), with n = (L - 1)/2, k = k w 


mod 1 = 1- modi (96) 

\ / U \ **V / a 

Thus if it was assumed incorrectly that the selected peak 
in the spectrum is associated with the set u , equation (96) 
will not contradict this assumption. It is therefore 
necessary that 


„ . „ / A* L-l Al\ 

+ ~2~T) (100) 

If the selected complex amplitude from the spectrum is 
incorrectly associated with set u, 

7ii^=7a Xjfl 


T mod 1 ?^0 

Equation (89) then yields 
k u mod 1 ;*0 


Using equations (99) and (100) (with u replacing u) in the 
preceding equation gives 


+N N K , n -(L~\)M/L 
«u- x a +is A/m+A l/L 


( 101 ) 


Note also that if T mod 1 t^O, in general the original value 
of (ojf/o3 r ) a mod 1 will also be changed. The only 
exceptions are when 


T mod 1 = 1 -2 



(=f)„ m ° dl< 5 


(97) 


in equation (94), and when 


where &j Kt n=ju Kn ~Ja Kn ' th e other hand (eq. (85)) 

C02) 

Equating now u Kifl and a K>ny given by equation (55), 
results in 



T mod 1=2 



0?)„ modl> 5 


(98) 


and 



it 



a 


in equation (95). 

Consider now the case when L is odd and 
n — (L-l)/2. Then the set u frequencies are associated 
with 


which is identical to equation (91). Hence, by comparing 
equations (101) and (89) and referring to equation (93), it 
follows that for the case n=(L- 1)/2 the previously 
derived equations, (94) to (98), will apply provided that 


N 

k = k u < ~2~ 



r - ^A, n -q- WTL 
Am + Al/L 


(103) 


and the set a with 
N (j3f 

J- 


Therefore, for n = n u = (L- l)/2 and x = /c a , 

= (F„ I 2 e -i2r[K a Av/LN+ (L - l)Al/2 L] = |g Q j /To,, 

where use was made of equation (76). Hence 

„ . „ / Av L - 1 Al\ 

ya', n + 2*Ja K .= ~ 2 *{ K a m + ~2~ T ) 


(99) 


This case is significant in particular for L = 1 (and 
therefore A / = 0) since it also applies for the discrete blade 
displacement field. 

Note that if the selected complex amplitude were 
incorrectly assigned to the set a rather than to the set u, 
subscripts a and u would need to be replaced by a and u in 
the preceding equations. These equations would then lead 
to similar results and conclusions as in the previous case. 

Having derived T for all n, it is necessary to obtain the 
range for Aj K so that appropriate r’s can be evaluated 
for any given choice of Am and A/. As discussed 
previously in connection with equations (92) to (98), it is 
desirable that T mod 1^0. Using equations (86), (92), 
and (93) (or alternatively (103), if n=(L- 1)/2), it 
follows that 
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-i<r/;v<i (104) 

This equation will be satisfied for appropriate T’s 
provided A j K n is within the following limits: 


- {Am - 1) <Aj K n <Am - 1 


- Am + A<Aj K>n -<Am- \ +A 


- {Am + 1) + A <j Krt <Am - 1 + A 
where 


A/ = 0 


A 


Al/L <1/2 
Al/L >l/2j 


(105) 


A = 0 n*{L-\)/2 


and T mod VO cannot be satisfied for all Aj K fJ . 
However, as will be shown later, even for L = 1 it is 
possible to choose transducer separation A3> so that 
T mod 1 5*0. 

Having determined conditions which guarantee that 
the selected amplitude will be correctly assigned to either 
set u or a , there remains to be determined the condition 
under which there exists only one j Kjn and therefore only 
one k in either set u or set a . 

Let (k u ) T and (a y/av) r in the following equations 
denote the true values for these variables. Depending on 
the values of j Uiin , there may be several other values of k 
and o)f/a) r satisfying the appropriate equations. Then for 
every such k = (k u ) f and coy/a> r = (oy/av) F , 


A — Al n=(L- l)/2 


(MK y n) T =^K >n ) F 


Note that if A/ 5 * 0 , IV 0 for the preceding Aj Ktn . 
However, if A/ = 0, equation (96) results for Aj K n = 0, and 
therefore no distinction is possible between sets u and a . 
(The distinction between these sets when A/ = 0 would be 
possible only if the upper limit in inequality ( 86 ) is 
replaced by N/2 .) 

Equation (93) or (103) should be evaluated for all 
possible Aj K n ’s given by equations (105). However, this 
procedure can be somewhat simplified at the expense of 
introducing another restriction. From equation (93) it 
follows that 

T mod 1 = ^ N(L Ay «." + A /) ] mod A " 

Av 

Therefore T mod 1 5*0 if 
[JV(Z, Ay K>/I + A/)]mod Av 

= [TV mod A v(L Ay^„ + A/)mod A^]mod Av 
= [TV mod A v(L Aj K n + A/)] mod Avr*0 (106) 

in view of equations (105). It is therefore necessary that 


or 

(^) r + ( * u)r +n u N = (^) + ( K u ) F +n u N 
so that 

-[(U) r -(“),] (l08) 

Because 0 <a if/w r <N, it follows from equation (108) 
that 

|(« H ) r -(*«) F |<JV (109) 

On the other hand, from equation (83) it follows that 

(«u) r - («u) F = \(ju K ,) T ~ (/«,.„) J = 6 J*- n Am + Al/L 

(HO) 

where use was made of the fact that ( y u ) = ( 7 ^ ) and 

* u K,n j* **K,n f 

the definition 


iVmod Av 5 * 0 

If in addition Av is a prime number, the product of any of 
the factors of N mod Av and (L Aj K>n + Al)<Av cannot 
equal Av> and therefore T mod 1 will not be zero. To 
summarize, T mod 5*0 if 

Al 5*0 


$/ K,n~ 0 u K r ) 0 u Kt „) , 


Consequently, in view of equation (109), 

^ — |<i 

Am + Al/L | 
and therefore 


TVmodA^O and Av = prime number (107) 

Note, however, that the last condition on Av is sufficient 
but not necessary. Equations (107) also hold for 
n= (L — 1)/2 when L is odd. If L = 1, Al must equal zero 


| 6 y K>w ]<Am-l A/ = 0 


IS/iml^Am A/>0 


(111) 


( 112 ) 
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Applying the mod 1 function to the absolute values of the ^ < y K>n < A- 1 

left and right sides of equation (1 10) results in ~ 2ir ~ N 


(118) 


|(O r -(O f | mod 1 = | J^% /L | mod 1 = 0 

because both k’s are integers. Thus multiple solutions for 
k will be ruled out if the right side of this equation is zero 
only when §/*,„ = 0. It is necessary therefore that 

< ii3) 

for the nonzero values of \t>j Kt „\ given by equation (112). 

To derive a simplified condition on Am and Al, 
analogous to condition (107), equation (113) is trans- 
formed as follows: 


in equation (83) or (82) results in 


N - 1 
N 



, A/ . N — 1 Av A l 

+ n L ~ Jl< ’ n ~ N L +n L 


(119) 


where again either subscript u or a could be used. If 
Al=0, equation (119) reduces to 


N - 1 N - 1 

- -jj- ( Am + 1) < -jj- Am (120) 

or, since Am <iV- 1 

- Am <Am - 1 (121) 


|(^n) 7* ( K u)/r mod 1 — 


(NL \&j Kl „ |)mod Av 

Av 


Therefore 

( NL | Sj K>n j)mod Av = [Mnod Av(l | dj K n j)mod Av J mod Av 

= [Amod A*(£|5y K>/I |)]mod Av = 0 

(114) 

where use was made of equation (111). Using the same 
reasoning as in the derivation of equation (107), it 
follows that if 

Nmod Avj*0 and Av = prime number (115) 

| Sj K n | must be zero and the solution for k within set u is 
unique. As for equations (107), the second requirement 
on Av in equations (115) is more restrictive than 
necessary. Note that equations (115) and (107) are 
identical except for the additional requirement, A/ 5*0, in 
equations (107). Conditions (114) and (115) also hold for 
set a and for /* = (L-l)/2 when L is odd. If L = 1 or 
A/=0, conditions (115) are replaced by 

N mod Am 5*0 and Am = prime number (116) 


For odd L, n=(L- 1)/2, and odd TV the limits on k from 
equations (29) to (32) are 
r 






mod 1 < ^ 


mod 1 > ^ 


( 122 ) 


(& 

(?).' 


. , 1 N ~ 3 

mod 1 < « , =- 

2 2 


, , 1 N - 1 

mod 1 > - , — 

2 2 


<K a <N-l 


(123) 


Therefore, using equations (99) and (100), for odd N 


N- \ ( Av A L-1A/^. 

“ ~W ( T + y * T 


1 N— 1 Av L-\ Al . 

— ± — — 7 “ 




2 N L 
mod 1 < 


2 L 




mod 1 > 2 » 


> 


(124) 


As mentioned previously, for the determination of k it 
is necessary to assume j. It is useful therefore to derive 
the plausible limits on j. Applying the limits on ay/a> r to 
equation (48) gives for both sets of frequencies 


1 N - 3 Av l L-l Al 

2 N L ± 2 L ~ Jk ’" 




- (N~ 1 (117) 

Where the positive sign is to be used for set u and the 
Substituting these limits and the following limit on y (for minus sign for set a. For n~ {L — 1)/2 and even N from 
either y - y u or y = y a ) equations (27) and (28) it follows that 
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2ir/N-b* 



and 


N-l 

N 



± 


L—1AI . 

2 L ~ jK ’ n 



Av L- 1 Al 
L ± 2 L 


(125) 


(126) 


&V 

(Note that A# and Am increase in the direction opposite 
to the rotation.) Hence one must introduce the corre- 
sponding phase angle correction in the transform of the 
second transducer. Therefore 

Is J e /f “‘= K)(®J*= (K 

xe n™ u ( Am + n/N)(a u y 


= (a u )(a- u ) e^rKv/N-S*) 


where the signs are to be used as in equation (124). If 
Al=0 and n= (L — 1)/2 


Am 


-A m<j K>n < < 


2 

Am - 1 


- 1 Am even 


(127) 


Am odd 


xe i2*-K u (Am + l)/N = q g'~*u K 


and 

+ (129) 

where y UK is the measured phase angle, which is distinct 
from the true cross-power phase y Ufc 


provided that Am <N/2 for even N and Am <N/ 3 for 
odd N. These limits are also applicable for the L = 1 case 
of blade displacements. 

The general procedure for finding k and oy/w r 
associated with the particular Q K rt can now be described. 
Given the measured values of y K>n , u K n (or a K n ) and 
(«// u) r ) mod 1, first compute the possible k u ’ s and K a *s 
from equation (83) or (82) using fs given by any one of 
equations (119), (121), (124), (126), or (127). Reject /c’s 
that are nonintegral and for each integral k compute 
u)f/a) r from equations (85) and (90). Reject those k* s for 
which (c Of/u> r ) mod 1 does not correspond to the meas- 
ured value or for which oy/w r does not satisfy inequality 
(86). The remaining k’s should all be associated with 
either the set u only or the set a only. If the transducer 
spacing was chosen so that T mod 1^0 for the 
appropriate T, there will remain only one k for which the 
associated a )f/o) r is between 0 and N. 

When sampling blade displacements and the trans- 
ducer separation A$ is such that 


The corresponding equation for the set a is 

cx) f \ /2tt _ Am 4- 1 

i).\» ) 2 ’‘‘— (130) 

The equations 


7 a. + 2 ir/, 


'a,- ^ 


a K> yu K ~~ ya K 

together with equations (14), (20), (129), and (130) can 
now be used to obtain 


\»r/ u WJ a Am + 8*/(2ic/N) 


(131) 


where Aj K = ju K ~Ja K - Comparison of equation (92) with 
the preceding then leads to the appropriate T, 


5 $ = A4> mod 



5*0 


(128) 


r _ »,A/k- 1 + b $/(2r/N) 
Am + S*/( 2 -k/N) 


(132) 


before taking the cross-power spectrum, one must com- 
pensate for the fact that at the time t = 0, when the blade 
m = 0 is aligned with the first transducer, no blade is 
aligned with the second transducer. The next available 
blade (i.e., m = Am + 1) will reach the second transducer 
after the time period 


For the same reasons discussed in connection with 
equations (92) to (98), 

T mod 1 j*0 

to correctly associate a particular complex amplitude 
with either set u or set a. 
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The limits on A j K in equation (132) can again be derived 
by using equation (104). Thus 

64> = 0 


- (Am - 2) < Aj K < Am 

- (Am - 1 ) < Aj K < Am 
- Am<Aj K -<Am 


54> 


1 

(2ic/N) ~ 2 


> (133) 


The possible nonzero values for dj K are given by 
inequalities (135). If 8$/(2ir/N) and A4?/(2tt/N) can be 
expressed by equation (134), it is sufficient to satisfy 
conditions (116) instead of (136). Equations (116) and 
(136) also apply for set a. 

In this case the estimates for j u and j a are given by 


<54> ^ 1 

(2t/N) > 2 


,N- 1.. „ . 1 /7V-1 \ .. 

" ~jt (Am+2) - J “'-2\'JT ) (Aw + 1) 


*\ 


If <5<l> = 0, for A/ k = 1, T = 0 and thus T mod 1=0. It is 
therefore desirable that 5$ not be zero (otherwise eq. (86) 
should be replaced by 0<coy/co r <7V/2). If 6$ is not equal 
to zero, equation (133) should be used to verify that 
T mod 1 is not zero for a particular choice of Am and 6$. 
However, if 5$/(2tt/N) can be expressed as a ratio of 
two integers, that is 


♦ 0 - 5 &) 


> (137) 


-( 


i-JL.) 

2ir/N ) 


N— 1 
N 


(Am + 2)<j a 


6<i> 


(2ir/N) 

and 

A$ 

(2ic/N) 


A[ 

L 


"\ 


•i(V) 


(Am -hi) 


Av 

T J 


(134) 


equations (107) will insure that T mod 1^0. (Note that 
here Al,Av , and L are used only to define 6$ and and 

that L does not indicate number of samples per blade 
passage during sampling.) 

As in the previous case, one must also insure that there 
is only one solution for k in a given set (u or a ). The 
expression 


These inequalities can be derived by using equations 
(122), (123), (129), and (130). 

The flutter frequency can also be obtained by using the 
individual blade spectra. The phase angle difference be- 
tween the two transducers based on the individual blade 
spectra of a single blade is related to the flutter frequency 
and the angular transducer separation as follows: 


•t>m-<t>m + 2ir/ 




*\ 


o^tm_ fo»ci y 


2tt 


(O r Uu) F Am+6$/(27r/A0 


V.-foJ ~(ju) (^) A * 


(138) 


J 


can be derived similarly to equation (110). To satisfy 
equation (86), 


where again the prime is used to associate a variable with 
the second transducer. The expression for the errone- 
ously assigned (to the set u) flutter frequency is then 


l3/> 


Am + 64>/(2x/N) 
and 

|6y K |<Am-l 

|S/J<Am 


<1 




N 


Am + 6$/ (2ir/N) 


Ay 


<5<i> = 0 


\ _[ 

\ A ' 1 

/ a L 

Am + S<&/(2ir/N) J 


N 


(135) 




where 5$ is given by equation (128) and Aj=j u —j a . 
Comparing equation (92) with the preceding gives 


Therefore k will be unique within a given set if A4> is such 
that 




Ay 


Am+8$/(2w/N) 


\N 


(139) 


r nHc 


Am + 8$/(2ic/N) 


mod 1 9^0 


8j K *0 (136) 


To be able to associate a selected complex amplitude with 
either set u or set a , it is necessary for 54> to be such that 
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r mod 1 ;*0. The ranges for A j to be used in the 
expression for Y can be obtained from equation (138) by 
applying inequality (86). They are 


1 < Ay <2 A/m — 1 
1 <Ay<2Am 
1 <Ay'<2Am + 1 


5<t> = 0 

d*/(2ir/N) <1/2 
5$/(2ir/N) >1/2. 


(140) 


Wu\n 

Am 4* 53>/(27r/A0 


jmod 1 




Since T = 0 when 5$ = 0 and Aj = Am, 5<l> should not be 
zero. Provided b$/(2w/N) and A4>/(27 t /N) can be 
expressed by equations (134), it is sufficient to satisfy 
conditions (107) to insure that T mod 1^0. However, 
because one now has to work with frequencies rather 
than nodal diameters, restrictions (97) and (98) apply. 

Assuming again the true and false values for ay/o> r 
(associated with the true and false values of j u ), equation 
(138) yields 


(144) 

From equations (138) and (86) the estimates for ranges 
of j u and j a are 

0</ tt <Am — 1 -Am</ fl < — 1 5<t> = 0 

> (145) 

0</ u <Am -(Am + l)</ ff < - 1 


(aA ~(^i) - bju 

\ 0) r/’f V^r/ f A$/2x 


Am+SQ/pir/N) 6ju (141) 


where 5y u = C/ w ) r -(/ w ) F - In view of the limits imposed on 

C0//C0 r , 



§iu 

Am + 5$/(2Tr/N) 


N<N 


(142) 


and 

|5/ u | <Am <54>>0 


\bj u \<Am - 1 54> = 0 


(143) 


Since the left side of equation (142) must be a positive 
integer, 

|(^) Lodi 

J \°>r/T y^r/F | 

= [ Am + ii/(2ir /N) ] m ° d 1 = ° 

Therefore, to have a single solution for oy/co r within set 
u 3 equation (136) must be satisfied with bj u replacing bj K . 
The nonzero values for |6y M | are given by inequalities 
(143). This procedure also holds for set a. If the angular 
separation between two transducers can be expressed by 
equations (134), condition (116) can be used instead of 
(137). As for T, for this sampling rate there is an 
additional restriction 


Table I (ref. 2) illustrates the phase relationship 
between the two pressure transducers 47.4° apart (a dis- 
tance that corresponded to Am = 5). The flutter point was 
the same as for figures 1 to 3, and therefore the true 
flutter frequency was 8.45 E. In the axial direction these 
transducers were located past the blade midchord posi- 
tion. In contrast to the leading-edge location, at this 
chordwise position there was a significant response for 
the n = 0 harmonic. Consequently a low-pass filter could 
be used to cut off the response beyond about ay/a> r = 25. 
This eliminated the higher harmonics, as well as the 
dominant part of the steady-state response and consider- 
ably simplified the analysis. 

In this case the true flutter frequency was known from 
the blade-mounted strain-gage outputs. Consequently the 
cross-power phase angles based on the known nodal 
diameters associated with each frequency in the first 
column could be computed and tabulated in the last 
column. Comparing these phase angles with the 
measured phase angles in the second column indicates a 
particularly satisfactory agreement for the backward- 
traveling waves. 

TABLE I. -PHASE RELATIONSHIPS 
BETWEEN TWO TRANSDUCERS 
47.4° APART 


Frequency, 

engine 

orders 

Cross-power- 
spectrum phase, 
deg 

Phase based 
on k, 
deg 

3.45 

122.0 

123.0 

4.45 

173.6 

170.5 

5.45 

215.8 

217.9 

6.43 

266.3 

265.3 

8.44 

256.7 

0 

11.46 

150.0 

142.1 

12.45 

201.7 

189.5 
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To illustrate the procedure presented in this section, 
flutter frequency will assumed to be unknown now, and 
the measured cross-power phase associated with the 
5.45 E response will be used to determine it. Note that 
A / = 0, Am = 5 is a prime number, and, for N~ 38, TV mod 
Am- 3^0. Therefore condition (116) is satisfied; 
however, condition (107) is not. From equation (121), 
-5<Aj K>n <4. For these f s equation (83) yields two 
values of k u that are close to being integers, 

k u = -3.04 3 (for j KtTl = -1) and k u = 34,96 — 35 (for 

j Kt n = 4). For these *’ s from equation (53) oy/w r = 8.45 
and -29.55. The last value associated with k u = 35 must 
be ruled out. From equation (82), the other two possible 
values of *’ s (corresponding to the same values of j Kt „) are 

*0 = 3.04-3 and k q = -34.96 35. For K a - -35 from 

equation (90) ay/ oo r - -2.45, and therefore this solution 
must also be ruled out. However, for k q = 3, 
oy/co r = 29.55, which is within the assumed limits given by 
inequality (86). Without further input therefore it is 
impossible to determine the correct solution. Had the 
determination of frequency been the primary object, a 
better choice for second transducer location would be the 
one that satisfies condition (107). In particular this 
condition requires that Ah*Q. 

Interaction with Steady-State Pressure 

The true fluctuating pressure on the blade surface must 
be defined in terms of the coordinate system (£,?j) fixed 
with respect to a blade, figure 5. The coordinate system 
( x,y ) in this figure is fixed with respect to the rotor, and 
the x coordinate is measured from the undeflected 
position of blade m — 0 along the outer circumference of 
the rotor. The instantaneous displacement of the m * h 
blade in the ( x,y ) coordinate system is given by bx m ; and 
the undeflected blade position, by x m . 

For the m th blade the pressure measured by a 
stationary sensor P t>m can be expressed as a sum of the 
unsteady pressure P m U-,t) and the steady-state pressure 
P that is, 



Figure 5. -Definition of (£,i?) and (x,y) coordinate systems. 


In the (x,y ) coordinate system, both P m and P s m depend 
on blade displacement 6x m . Expanding equation (148) in 
powers of dx m to first-order results in 


Pt,m x m>t) +P s,m( x x m ) 


9Pm 

dx 


6x~0 


dP I 

Xdx '”- ± J? L \ bX m ( 149 ) 

ax \8x=0 

\x—x m \< <Riz/N 

where the derivatives are evaluated for dx m = 0. Since it is 
expected that P m < <P s ,m> 

dPm ^ ^ dP s m I 

bx 8x = 0 dx ISx = 0 

and 

dP I 

Pt,m~Pm( x ~ x m>0 P St m( x ~~ x m) I $ x m 

\8x = 0 


P t ,m=Pm(U)+Ps,m(S) |£|<<^*^ (146) 

where the condition on £ implies that the pressure distri- 
butions are restricted to the immediate neighborhood of a 
blade. Both positive and negative £’s are allowed. It is 
also assumed that the blade displacements are much 
smaller than the blade spacing. 

After the coordinate transformation 

*-*m = £ + 5 *m(0 (147) 

is introduced, equation (146) becomes 

m( x ~~ x m — ^ x m> t ) + P s>m (x—X m — & x m) (l 4 ^) 


(150) 

Transformation into a stationary reference frame by 
using Jhe real part of equation (11) with ^k—^bx,k and 
and the real part of a replaced by bx m results 
in 

dP \ N ~ l 

P'-PM+PAri-^l E 

aX \8x=0k = 0 

xcos +Ar ) 77 (m+ ^ Ar) + (151) 

\p- ti m \<:e< < or \n~ mL\<e+qLN< < +gLN 
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where is defined by equations (36) and (37). When 
L= 1, P(ii) is given by the real part of equation (11) (in 
which case P, P s , and ( dP s /dx ) | 5x= o become functions of 
m + qN ). When L> 1, P(fi) is given by the real part of 
equation (47). The third term in equation (151) is 
proportional to the steady-state gradient evaluated for 
5x m = 0. In the ( x,y ) coordinate system, it represents an 
apparent unsteady pressure component that is created by 
the displacement of the steady-state pressure field. 
Because the steady-state pressure gradient is expected to 
be much smaller away from the immediate neighborhood 
of a blade, the contribution of this term can be neglected 
to a first approximation in this region. Therefore 

P t ~P(li) +Py(/*) for |/* — mL\>e + qLN, L>\ 

Also, to properly use equation (151), it is assumed that m 
is incremented in the last term when m = ~ I e I and 
not 0. The pressure P is now defined for all ranges of fi 
when L > 1 . 

If data are sampled over Q revolutions and then 
arranged into a matrix so that each row starts at a fixed 
v = vq (£>1) or m=mo (£ = 1) and contains either LN 
(L > 1) or N (L = 1) samples, the phase angles associated 
with each row (i.e., each revolution) will be random for 
the first and third terms in equation (151) and constant 
for the second term in this equation. This follows from 
the fact that all harmonics for the first term (eqs. (1 1) and 
(47)) and the third term contain a nonintegral shift 
±uf/o) n whereas the second term originates from the 
steady-state term and therefore all harmonics for this 
term are of integral engine orders. Therefore, by 
averaging Q rows, column by column, the contribution 
of first and third terms will be essentially eliminated and 


only the contribution of the second term in equation 
(151) will remain. This presents a way of eliminating the 
steady-state contribution from the spectrum. Averaging 
is most conveniently performed in the time domain by 
using the once-per-revolution trigger to mark each 
revolution (i.e., fix v 0 or mo). (Note that revolutions for 
this procedure do not have to be consecutive.) Once the 
steady-state distribution is obtained in the time domain, 
it can be subtracted from the total signal for each q and v 
(or m) leaving only the unsteady components in the 
signal. The same procedure is applied for both pressures 
and displacements with the exception that the third term 
in equation (151) is nonexistent and only the L = 1 case is 
possible for displacement signals. 

For the pressure signals it is also desirable to eliminate 
the contribution of the third term. To do this according 
to equation (151) would be tedious. However, as will be 
shown, its presence does not affect stability criteria, and 
consequently it is not an important contaminant in the 
spectrum. Also, for £>1, when it is desired to obtain 
pressure amplitude distribution in the blade passage, L is 
usually selected sufficiently large so that the proportion 
of samples falling in the range |/x - fi m | < e can be counted 
to be small. Some other ways of minimizing the contri- 
bution of this term when blade surface pressures are 
desired and L= 1 are discussed in reference 3. 

Figure 6 illustrates a linear pressure spectrum for which 
the steady-state response was not removed before spectral 
analysis. The flutter point is the same as for figure 4. 
Comparing magnitudes of the unsteady and the steady- 
state amplitudes justifies the assumption made in de- 
riving equation (150). 

Some experimental evidence regarding the contribution 
to the unsteady pressure due to the displacement of the 
steady-state pressure field can be found in reference 3. 



Frequency, engine orders 

Figure 6. - Linear spectrum of leading-edge pressure transducer; integral engine-order harmonics retained; engine B. 
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( 156 ) 


Computation of Aerodynamic Work 

In case flutter is mistuned and several modes are 
present, it is of interest to establish which are the least 
stable modes. If both pressures and displacements are 
measured, it is possible to do this by computing the aero- 
dynamic work for each mode. For completeness there- 
fore the two-dimensional expressions for aerodynamic 
work will be derived by using the modal expansions for 
pressure and displacement (as would be obtained from 
the overall stationary spectra for these variables). The 
most frequent case of torsional flutter will be assumed. 

Taking the positive direction for the moment to coin- 
cide with the positive sense of variation of the 0 
coordinate in figure 7, the contribution to the moment 
(per unit span) from the t/f element of the chord for the 
m th blade is 

dM m = AP m ( X c-n rff (152) 

where A P m denotes suction- minus pressure-side 
unsteady differential pressure, \c is the distance from the 
leading edge to the torsional axis, and the f coordinate is 
measured from the leading edge along the chord. The dif- 
ferential work can now be obtained by multiplying equa- 
tion (152) by the differential displacement 

dO m = \d(6 m + 6- m ) (153) 

which results in the differential work 

cVW m = I AP m d ( 0 m + fr m )( X c - f) rff 

= ^ A P m e io3 f ! \6 m | (e‘ ( u /‘ + *0 - e ~ '(“/* + 

x(xc-f) d&(u f t) 

= ^AP OT pT w |e-^ m (e 2 'V+^ m )-l) ( X c-n d(a f t) 

= ^AP M C(e 2 'V+^J-l)( x c-f) d&lujt) (154) 

where 4>Q m is the phase associated with 6 m . Integrating this 
expression with respect to j6 = ay/ + cf)$ m from 0 to 2ir gives 

dw m = La PM xe-n df j*' d& 


N-\ 

dW=-iw( X c-t) off D A P m e' m 

m = 0 

The summation in this equation can be expressed in 
terms of DFT’s of A P m and 6* m > which are denoted by 
A &p y k and Afik. Using equation (12) in terms of these 
variables^ with k substituted for the dummy index k in the 
sum for d* m , it follows that 

N - 1 _ N-\ N~\ N-\ 

D L A APik Al K e^m(k- K )/N 
m ~ 0 m = 0 A: = 0 k = 0 

N— 1 N— 1 _ N— 1 

= £ E A AP,kA'e,K E e‘2Tm{k-K)/N 

k = 0 k ~ 0 m-0 

(157) 

Since 

n - i r=o k^K 

F eilTtm ( k — K ) /N s 
m = 0 U N k = K 

equation (157) becomes 

N - 1 _ N— 1 _ 

E A P m e' m =N E A AP,kAe,k (158) 

m = 0 k=0 

Substitution of this identity in equation (156) results in 

N~ 1 _ 

dW=-iir( X c-n d£N E A AP,k A 6,k 

k = 0 
N~\ 

= -/tKxc-O d£N E 'l^AP.kl \A e , k \ e 1 ^-^ 

k~0 

(159) 



Figure 7. - Location of measurement ports. 


= -i*APMxc-n df (155) 

Summing equation (155) over m results in 
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and 


(Q/),= 


l 


N- 1 


dfV=Ti(xc—f)N ^2 Map, A: I sin(^/>*— 0^*) 

k = 0 


(160) 


Integration with respect to f gives 

ir=w7v r (xc-n Map.atI M<ul 

J0 |_£ = 0 


x sin(^Ap 


'.k-'Pe.k) J 


d$ (161) 


or for any particular k 

W k = irN f (xc-DMap,*! s > n (^AP,/t — 

•1 0 

(162) 


From equations (156) and (159) it jfohows that only the 
imaginary parts of the products A P m d* m and 
contribute to aerodynamic work. Frequently the phase of 
pressure (or moment) is measured relative to the displace- 
ment. Then one can state that only the imaginary part of 
pressure (or moment) contributes to the work. Returning 
now to equation (151), it can be seen that the third term 
in this equation is either in phase or tt radians out of 
phase with the displacement (i.e., it is real). It follows 
therefore that this term has no significance for computa- 
tion of aerodynamic work and is irrelevant for stability of 
the system. 

Equation (162) is frequently stated in terms of dimen- 
sionless imaginary moment coefficient r ) / ; that is, 

w k= c2 ^\A e>k \ 2 (C M>k ) I (163) 


where p is the density, V is the relative velocity, and 
( C\f,k . ), is given by 


( c M,k),~ 



A AP,k 

P V 2/2 


x sin ('PbP'ic—'l'e.k) 



(164) 


The logical extensions for a mistuned rotor are 


Map.xI 


N- 1 

*5o 


xsin (.'pAP, k -h,k) 


L 1 \Ae,k\ 2 ( c M,k), (166) 

N-l *=° 

£ M 

k = 0 


Concluding Remarks 

Analytic expressions were developed that describe the 
displacement and the pressure spectra in a stationary 
reference frame. The sampling rates considered were the 
once per revolution, the once per blade, and for the 
pressure, a multiple of the once per blade. The once-per- 
revolution sampling involved multiple aliasing and hence 
appears to be less desirable than the once-per-blade 
sampling, which resulted in a spectrum that is folded only 
once. Despite the folding, all A possible modal responses 
can be easily obtained from this spectrum if the true 
flutter frequency is known. 

Analytic expressions were found to be essential for 
interpreting a complete pressure-field spectrum, which 
results when sampling at a rate multiple of the blade 
passing frequency. A procedure was described that allows 
the unsteady pressure field to be transformed into a 
rotating frame of reference yielding the distribution of 
the unsteady amplitude and phase around the rotor. For 
this spectrum, also included were the expressions for the 
nodal diameters. Although the complete pressure-field 
spectrum was found to be somewhat more difficult to 
analyze, its advantage is that there is no need to locate the 
blades in the sampled time domain record. (This step is 
necessary for the once-per-blade sampling of blade sur- 
face pressures.) 

Common to all stationary spectra is the fact that the 
true flutter frequency cannot be determined with a single 
sensor. A procedure was described that allows its deter- 
mination by using two sensors. All three sampling rates 
were considered and the aliased parts of the spectra were 
included in the analysis. 

Equations were derived that illustrate the contribution 
of the steady-state field to a stationary pressure 
spectrum. 

By using the expression for aerodynamic work, it was 
shown that stability of the system can be described in 
terms of the stationary cross-power spectrum between the 
differential pressure and the displacement. 


n V2 JV— 1 

W= p —-c 2 *N(C m ) ! £ \Ae, k \ 2 
* k - 0 


(165) 
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Appendix - Interpolation Formulas 


The formulas presented in this appendix were 
originally derived by Donald C. Braun of Lewis. They are 
presented here for completeness. The object of these 
formulas is to relate the amplitude A g , the phase Og, and 
the frequency f g of the cosine function 


Assume first the rectangular window w(t)= 1; then 

1 

' ± 1 +/? 

(A2) 


a(t) =A g cos(2ir/gt + 0g) (Al) 

to the parameters obtainable from the nearest frequency 
lines in the measured spectrum. (Note that the symbols in 
the appendix are defined independently of the main text 
and the symbol list.) However, because different window 
weighting functions are frequently used in the spectral 
analysis, the preceding cosine function will be modified 
as follows: 

h(t) =a (t)w(t) =A g cos(2 Trf g t + d g )w(t) 


s sin(xu) 1 1 

and 


T(Af-l) 

°g~Vr v 


(A3) 


Thus equations (A2) and (A3) give the A g and 6 g required 
in equation (Al). 

Consider now the Hanning window given by 


where w(f) is the weighting function. The DFT of h{t) is . r . _A h / 2irn\ n 

defined as w{nis) ~ T V ~N ) 


i N-l 

H k =nrr £ h{nT s ) e-^knlN 

iV n = 0 


where A h =y/ 8/3 is used to normalize the mean square 
value of the window. Then 


k = 0 N- 1; n = 0 N - 1 


v~ ± 


2-R 

l+R 


where T s is the sampling interval. 

If f r denotes the frequency of the dominant frequency 
line (i.e., closest to the true frequency f g ), and /3 the 
frequency bin spacing (/3 = l/(NT s )), the interpolation 
fraction v is defined by 


v = 



1 1 

- - <u< ^ 
2 2 


xr 2(l - u 2 )sin(27rt;/A r ) 
g Ah sin(w) 


\H r 


and 


(A4) 


O g ~<t> r —irv (A5) 

where in the expression for A g the approximation 


The amplitude ratio is defined by 



where the positive sign is used when_|// r+ i|> 
and the negative sign when |7J r _i[> |// r +i|. 


tan(7 w/N) 

V tan {t/N) 

was made. Equations (A4) and (A5) give the parameters 
required in equation (Al). 
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